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THERMODYNAMIC ANALOGIES FOR A SIMPLE 
DYNAMICAL SYSTEM 

By Edwin Bidwell Wilson 

Introduction. In a recent article * I have set forth some of the most 
elementary notions in regard to the statistical mechanics of dynamical systems. 
The object of the present communication, which may be regarded as a sequel 
to the other, is to investigate a simple dynamical system for the purpose of 
exhibiting certain analogies between its statistical properties and the elemen- 
tary thermodynamics of gases. The dynamical system which I shall choose is 
the somewhat familiar household article known as the Japanese or Chinese 
dinner-gong and consisting of a series of gongs of varying sizes suspended 
one above another upon a vertical string. The number of degrees of freedom 
of a set of several gongs is considerable ; and to simplify the work it will be 
assumed that the only motion which the gongs may have is one of rotation 
about the vertical string as axis. Under this hypothesis the number of degrees 
of freedom will be equal to the number of gongs in the set. 

To find the equations of motion of the system in the Lagrangian or 
Hamiltonian form, it is first necessary to write (}own the kinetic and potential 
energies of the system. If I^, I^, I3, • • • are the moments of inertia of the 
successive gongs beginning with the top one and if qi, q^, qs, • • • are the 
angles through which the respective gongs have turned from the position of 
equilibrium in which there is no twist in any of the strings, the angular veloc- 
ities of the gongs will be respectively qi, q^, q-i, • ■ • and the kinetic energy 
of the system will be 

T=i{Ixq{+I,ql + I,qi+---). 
The momenta corresponding to the different coordinates are 

p,= - = I,q, or q, = j^; 

• Applications of Probability to Mechanics, Annals of Mathematics, volume 10, pp. 
129-148 (1909). References to this article will be given simply by the page numbers in pa- 
rentheses without the quotation of the title. 

(149) 
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and hence as a function of the momenta the kinetic energy is 

For the purpose of getting rid of the constant w in some of the subsequent 
averages, the assumption will be made that the product 

and the kinetic energy will be taken in the form 

T=€p = -rr(aipl + aiPl + aapl+ ■ • •), a^a^as ■ ■ ■ = 1. (1) 

The potential energy of the system will consist of the energies stored in 
the twisted strings. The twists in the successive strings beginning with the 
top one and running down are 

<?i. ?2 — ?n ?8 — Ja* • • •• 

The potential energy will be proportional to the square of the twist, where 
the factor of proportionality may depend on a variety of things such as the 
kind of string, the diameter of the string, the length of the string, and the 
tension in the string due to the weight of the gongs. In order that there may 
be some possible variation of the system as a whole it will be assumed that the 
length of the uppermost string is adjustable ; this could easily be arranged by 
running the string between two closely pressed wheels, which could be moved 
up and down along the string, instead of merely pinching it in a vice. It is 
then important to know at least how the factor of proportionality varies with 
the length of a uniform string when nothing but the length changes. Now- 
from the symmetry of the figure it appears that to twist a string through an 
angle q is equivalent to twisting each half of it through the angle i q in oppo- 
site directions. Hence if the factor of proportionality is supposed to vary 
with the 7ith power of the length and is written as kl", the equation of the 
whole potential energy to the energy of the two halves gives 

^/»,-2^@"(f)^ or l = i(i)" 

* The use of ty for kinetic energy, «? for potential energy, and« = T + V= ep + e,loT total 
energy is adopted to conform to the usages of Gibbs's Elementary Principles in Statistical 
Mechanics. 
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and it is seen that n = — 1. Tlie potential energy of the system will there- 
fore be written in the form 

F= e, = kJlYi + *j^i"'(?2 - ?i)' + hlj'iqs - ?2)' + • • •• 

For the sake of simplification let it be assumed that an adjustment between 
the sort of string and the diameter of the strings and the units of measure- 
ment has been made so that 

kfl^'^ = v, k^lj^= IT, ■ . . ; and let kili^ = ■jrv~*, 

so that the parameter v is proportional to the fourth root of the length of the 
uppermost string of variable length. Then 

F= e, = ^lv-*ql + (q^ - q^y + (q, - q,y + • • •] 

is the expression for the potential energy of the system. 
The diiferential equations for the motion would be 

Pi +(v-* + l)qi -q^ = 0, p^+2q^-qi-q3 = 0, • . • 
or 

Pi+(^*+l)?i-?i! = 0, gi = aiPi, Pi+2qi-qi-q3=i0, q-^^a^p^, ■ ■ • 

by direct substitution in the formulas of the Lagrangian or Hamiltonian equa- 
tions (see page 147), where 

T- V= ep-€g = -rrlaipl + a^pl + a^pl + v-*q\-(q^- q^y 

or 

T + F= €p + e, = 7r[ai/i + a^pl ■{■ ci^pl -\. ■ ■ ■ + v-*' q\ + (q^ - qiY, 

If D denotes differentiation with respect to the time, the Hamiltonian equa- 
tions after the elimination of the momenta are 

[li)^ + (ir-*+l)]?,-?, =0, 
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Hence, if for example the number of gongs is three, the solution of the dif- 
ferential equations* depends on the solution of 



-Z>2 + v-*+l - 1 



-1 





-Z>* + 2 
- 1 




- 1 



-D^+ 1 



= 0, 



which is an algebraic equation of the third degree in IP. If there had been 
five gongs the solution would have depended on the solution of an algebraic 
equation of the fifth degree in D^. As these equations are necessarily literal 
when the length of the uppermost string is adjustable and v cannot be assigned 
a numerical value, the solution of the algebraic equations is impossible in 
practice. This will make no difference to us inasmuch as the discussion of 
statistical mechanics is based on the differential equations and not upon their 
integrals. 

1. The canonical distribution. Suppose for the moment that the 
system consists of three gongs. The state of the system is at any instant 
completely determined by the values of the six quantities g-j, g^, q^, pi, p^, p^. 
The representative point of the motion would therefore be located in a space 
of six dimensions. Consider an ensemble of these dinner-gong systems dis- 
tributed with the density or better with the probability (see page 148) 



^_. 



P = e ® , € = 7"+ F= Cp + €y, -^ and constant. 



where 



€ = TT [a^pl + a^pl + a^pl + i' *?! + (S'j - ffi)* + (^s - ^a)*]- 

The distribution will be in statistical equilibrium because the probability, or 
density which is proportional to it, is a function of the energy alone and the 
energy is an invariant of the motion. Let it further be assumed that the dis" 

* See Forsyth's Differential Equations, 3rd edition, p. 805. Or Cohen's Differential 
Equations, p. 150. Or any treatment of simultaneous linear differential equations with constant 
coefficients. 
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tribution covers the whole of the six-dimensional representative space so that 
the relation 

— OO 

must be satisfied (see pages 132, 142). This distribution will be considered 

as the canonical distribution and averages over it will be investigated.* 

The equation (2) may be used to determine the constant ^ in terms of 

the other constants of the problem. For as ylr and @ are constants, the 

* 
factor e* may be passed out from under the sign of integration and may be 

transferred to the other side of the equation. Then 

— 00 

e <H® dpi / e "'^ cZjOg / e "^e dp^ X 

• « J — 00 J — 00 

. « '^dq,J_ e- « d,,J_^ e- e dq,. 

Here the six successive integrations are to be carried out by beginning with 
the one for q^ and working backward through the list. As gj is regarded as 



• As a general comment on the assamptlons underlying the foregoing use of a canonical 
distribution, one point In particular should be mentioned. It has been assumed that the repre. 
sentatlve points for the motion of the systems In the ensemble are spread oyer the whole of the 
representative space, that Is, that the coordinates q^, q,, q^, . .. and momentup,, p^, p^, . . . may 
have any value. This necessitates the admission that the angular velocities of the discs and the 
twists in the strings may be indefinitely great. It goes without saying that this is physically 
Impossible : for If the strings were allowed to twist very much, the law of torsion would de- 
part from the simple linear (Hooke's) law and the expression of the potential energy would 
change ; furthermore the twisting would perhaps materially shorten the strings and Introduce 
complications on that score. It should be noted, however, that as the coordinates and momenta 
Increase, the energy increases with their square, and the probability assumed for the distribu- 
tion decreases exponentially with the energy. Hence the probabilty of systems with large ve- 
locities and twists is almost negligibly small compared with that of systems with small velocities 
and twists. The errors made by entirely neglecting tlie systems In which the velocities or twists 
were so large as to introduce physical dlfflcultles would therefore probal)ly be small and could 
certainly be made small by restricting 9 to sufficiently small values. Hence there Is no need 
of worrying very much over such difficulties. 
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constant for the integration with respect to q^, we may write 

/+0O ir(g3 — gj)' /•+«» T(ga— gj)' \J% I— -— 

e » <i?3=/ e « d(g'3-g'2) = -y V/7r= V0- 
•00 J— 00 V"" 

The result is independent of q^ ; in fact the integrals may be regarded as 
multiplied instead of iterated. The result is 

the last reduction being due to the assumption (1 ) . Hence 

Vr = - 2@ logv - 3e logo. (3) 

It is readily seen that if there had been n gongs, the expression for i^ would 
have been 

A|r = _ 2 © log V - n © log @. (3') 

Another important quantity is the index t) of probability which is defined 
as the logarithm of the probability. Thus 

^ = logP=^!^, P=e\ 

The average value of the index of probability jj taken over the ensemble is 
found by evaluating the integral 



r,= I • • . ve^'dpi ■ ■ . dqs. 



This may be evaluated directly ; but the work may be very much shortened 
by adopting the device of differentiating (2) with respect to ® regarded as 
a parameter. Upon the left-hand side, the differentiation may by Leibniz's 
rule be carried under the sign of integration. The result is 

Now @ may be canceled out and the result transposed so that 

a@Xv ■]'' ^^^■■- ^^3=j_-; -j-^e^ dp,...aq, 
or _ . 1 =jf^.^. .J r,e^dpi ■ . . dqs = V 
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That the factor d->^i'd @ may be transferred to outside the integral sign 
depends on the fact that ^ in (3) does not contain the variables with regard to 
which the integration is performed. 

By differentiating i/r with respect to @ the value of ^ is 

^ = - 2 log v - 3 log e - a. (4) 

And now the average value i of the energy e may be found very simply. For 
as 1^ and @ are constants, it follows from 

V = -@- that ,, _ -^- . 
Hence i = ilr - 6 ^ = 36. (5) 

Moreover if there were n gonga, the results would be 

77 = — 2 log V — n log @ — n, (4') 

e = -s^-%-ij^nB. (5') 

The parameter v which enters into the expression for the potential energy 
and which determines the length of the uppermost string may be regarded as 
an external coordinate of the system. The difference between an external 
and an internal coordinate is simply that the internal coordinates are those 
which take an active part in the motion, and ih the solution are to be ex- 
pressed in terms of the time ; whereas the external coordinates fix the system 
as a whole and take no active part in the motion, although different motions 
will arise from varying such coordinates. Now the forces are derivable from 
the potential energy and they will be of two sorts, internal and external. If 
Qi, Q2, Q3, • • • are the internal forces and A is the external force,* then 

* The external force A should not be confused with the force which is required to support 
the system of gongs from falling under the influence of gravity. The force of gravity does 
not enter essentially into our problem. In fact it might be convenient to imagine that the force 
of gravity were continuously diminished toward zero as a limit as if the system were removed 
to a point indeflnitely near to center of the earth. It may be noted, however, that the force of 
gravity does no work on the system even when the length of the uppermost string is changed : 
for the assumption has been made that the wheels which pinch the string are moved up and 
down to regulate the length rather than that the string should be paid out between them 
In sach a way as to raise or lower the system. 
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It is easy to see without recourse to formulas that if the uppermost string is 
lengthened, the system loses energy : for a given angular twist upon a short 
string represents a greater amount of energy stored in the string than does the 
same twist in a longer string. 

The average value of the external force A over the ensemble may readily 
be computed. It is merely 



_ f+m r a. >!' — ' 



To evaluate this the device of differentiating (2) will again be used, but in 
this case the differentiation will be with respect to v. Note that yjr contains v but 
@ does not, and that as Cp does not contain v, the derivatives of e and e, with 
respect to v are equal. Then 



ix-r^ 



ci-fe-y"^*.-"^'-"' 



Hence ^ = -^.1 = 2-. (6) 

dv V . 

The result is independent of the number of gongs. 

2. A tlierniodyTiainic analogy. Let^, F, T, U, aS be respectively 
the pressure, volume, temperature, energy, and entropy of a gas. Let the gas 
be assumed to be in so far perfect or ideal as to satisfy Boyle's characteristic 
equation p V— iJT without sensible error. Further let the units of measure- 
ment be those ordinarily adopted in theoretical and practical physical chemis- 
try. • Then jj F = 2 T is the form of the characteristic equation. Moreover the 
energy of a perfect gas varies directly with the absolute temperature and is 
C/"=c„r where c„ denotes the specific heat at constant volume. In the 

♦ At the absolnte temperature T = 273° centigrade of melting ice and under tlie normal 
atmospheric pressure of 76 cm. of mercury or 1033 gr. per sq. cm., the gram-molecular volume 
is 22,380 cu. cm. Hence 

^ = '273 X ISao = ' (or Just a trifle less). 

See Walker's Introduction to Physical Chemistry, 2nd edition, p. 30, or any similar work. 
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physico-chemical system of units c^ = 3 for a monatomic gas ; for the ordinary 
diatomic gas c^ = 5 ; and for more complicated gases the values are 6, 7, 8 or 
numbers in the immediate vicinity of these integers.* Hence there is the 
equation 

U=c„T=nT, n = 3, 5, 6, 7, 8, 

for the energy of any ordinary simple perfect gas. The entropy and free 
energy of the gas have the familiar forms 

S = c,log T + 2 log F + const. and ^ = - 2 Tlog F - c^Tlog T; 

and finally, the thermodynamic surface for the gas, which is the equation 
between the energy, entropy and volume, is 

V' F^ = Ee^y Kaa arbitrary const. 

With these data for gases in mind it is easy to set up an analogy between 
the thermodynamics of a perfect gas and the averages obtained over a canoni- 
cally distributed ensemble of vibrating dinner-gongs. The results may be 
exhibited most graphically in double columns. 

Tkermodynamica of a Statistical Mechanics of the 

Perfect Gas Dinner-gong 

temperature T constant (modulus) © 

energy V average energy i 

negative avemge index 
entropy 8 ^^ probability " ^ 

free energy i/r constant i/r 

volume F external coordinate v « Z* 

average external generalized 
pressure i, = - (11)^ force ^= - (||:)^ 

specific heat c, degrees of freedom n 

equation pV= 2Tt equation Av = 20 

* These nambers are of coarse the "degrees of freedom" of the molecule. Chlorine 
is a somewhat anomoloas diatomic gas with c, = 6. Ethylene give c„ = 8. 

t This is the so-called characteristic equation of Boyle for perfect gases. 
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equation U= c^T equation i = n@ 

equation •^ = — 2 Tlog V — c^Tlog T equation •^ = — 20 log v — n@ log @ 

equation S = c„log T + 2 log F equation — rj = 7ilog@ + 21og« + n 

+ const.* 
equation V^ V^ = J£e*t equation e^v* = k e^^ 

It therefore appears that the statistical mechanics of the system contain- 
ing n gongs has all the formal properties of the thermodynamics of a perfect 
gas except that the arbitrary constant which enters into the entropy of a gas 
(and consequently into certain other formulas) is determined in the dynamical 
system. It may be well to note in particular some of the correspondences. 
The volume of the gas is related with the fourth root of the length of the up- 
permost string. If the volume is increased, the gas does work and loses 
energy ; if the string is lengthened energy is lost to the system. To the 
specific heat corresponds the number of gongs, that is, the number of degrees 
of freedom in the system — and this is precisely the number of degrees of 
freedom assigned to the molecules of the gas. As the above equations 
contain the basis of all the thermodynamics of the gas and have identical 
analogs in the averages taken over the canonically distributed ensemble of 
dynamical systems, it appears that the development of the two could never 
diverge. This conclusion, however, is a bit hasty for the reason that nothing 
has been developed which is analogous to the mixing of two gases. As a matter 
of fact the mixing of gases depends on a principle in addition to those here 
given for a simple gas and contains a somewhat troublesome paradox known 
as Gibbs's paradox. The matter of mixing gases and the analogous combina- 
tions of dynamical systems will therefore be omitted from the present 
discussion. 

There is one point of general import which should be mentioned. The 
reason for replacing the study of one particular vibrating system of gongs by 
the study of average values over an ensemble was that nothing was supposed 
to be known concerning the initial conditions of any one system (page 

* The entropy necessarily contains an additive constant as It is defined by an integral for 
which there is no natural lower limit — the absolate zero of temperature being escluded because 
the integral diverges for T = 0. Similarly ^= U — T'iS'contains an arbitrary constant and may 
be written ^ = — '2T log V — c» Dog T + (c» — const.) T; but the origin of entropy may be so 
choseu as to mal<e this final term vanish. The constant n which occurs in — ^ Is not arbitrary. 

t The equation of the thermodynamic surface. The constant K is arbitrary owing to 
he additive constant In S; but k = «— ". 
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137-8) ; the only known fact was the total energy of the system and the 
differential equations of motion immediately derivable therefrom. The ques- 
tion naturally arises : How nearly do the average values represent the condi- 
tions likely to obtain in a system selected at random? This is quite like 
the analogous question in the theory of errors of observation : How re- 
liable is the result as determined by the theory? A measure of the 
reliability of the result is assumed to be found in the value of the mean square 
error ; and for the dynamical systems the average value of the square of the 
departure e — i of the energy from its mean value may naturally be taken as 
a measure of how far the average state differs from the state of a system 
selected at random. It remains, therefore to evaluate the average of (e — e)*. 
Now it is clear that the relations 



(€ - i)* = 6* - 2€€ + ? = e* - 2ii + i* = e^ -i^ 

hold for the averages. It is therefore necessary to evaluate ? ; and the device of 
differentiating under the sign of integration may again be advantageously em- 
ployed. By definition and by (5) we have 

• • • / ee ® dpi • • • dq3 = 30. 
Now differentiate with respect to on both sides. Then 

• • • / €*e ® dpi • • • dqi-yjr . . . ee ^ dpi ■ • • dq.^ 

+ @ gljfjj -J^e « dp,. . . dq, = 3&K 

or ?_^i + @Ui= 3es. 

The substitution from (5) and (3) or from (5') or (3') gives 

? = 12@'^ or n(n + 1)0* ; ? - ? = 30^ or h02, 
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according as the number of gongs is 3 orn. Hence the ratio 



(e-i)*: i*= 1:3 or l:n. 

It therefore appears that as far as the energy is concerned, the mean square error 
in assuming that the average represents a system selected at random is consider- 
able when n is small but not very large when n is great and becomes quite 
negligible when n is very great. A similar calculation could be carried out tor 
the quantities 

(t) - vy and (A - Ay ; 

but this will be left to the reader. 

A simple device will give a system in which the average condition better 
represents the condition probably found in a system selected at random. For 
suppose that instead of considering single sets of n gongs, the attention is 
turned toward pairs of exactly similar sets or in general to a system consist- 
ing of iV of the single systems of gongs. This compound system will have Nn 
degrees of freedom if each individual set has n gongs. The energy of the 
compound system will be 

1 1 

[a^pii + Pw +pli + • • ■ + v-*^u + {qii - quY +i9si - qnY + •••]• 

If the compound systems are to be canonically distributed over their repre- 

sentative space of 2Jfn dimensions, the probability will be e * where « has 
this new value. As there are no relations between the diflferent sets of coor- 
dinates and momenta pu, Pa, Pah • • •« Qw 9ih Qzh • • •. i* is readily seen that 
the quantities connected with the compound systems are 

-f = — 2iV ®\ogv — 27n & log © 
rj = ^ = — 2i\riog V — Ml log — iVn 

e = yfr — &ij =i J\7w0 
8i^ ^ ,^ 



(e _ i)2 : ? = 1 : iV^n 
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By taking iV suflSciently large, the probability of the energy of an arbitrarily 
selected system being noticeably at variance with the average energy is very 
small. Moreover the only difference between the new and the old equations is 
that ■^j Tj, €, A, are multiplied by the constant iV" — a difference which is of 
the same sort as if these quantities were measured in different units. The 
thermodynamic analogy therefore persists in all its essentials. It even has an 
additional advantage in that now any one system represents the average or vice 
versa with great probable precision, which is surely true of any particular 
sample of gas and the average of all similar samples. 

3. Microcanonical ensembles. A microcanonical distribution is 
derived from the canonical by confining the attention to those systems of the 
ensemble which are characterized by representative points between the surfaces 
e = C and e + de = O + dO in the representative space. Now as the proba- 
bility depends only on the energy and the energy is sensibly constant between 
these two surfaces, it follows that the actual infinitesimal probability (see page 
134) that the representative point shall lie between the two successive 

surfaces is 

tizl ^^ t:=J/dv\ 

where V represents the volume of the six-dimensional space included by the 
surface e = C — a quadric surface in pi, p^, p^, qi, q^, q^, — and d Frepresents 
the increment of volume between this surface and the succeeding one. 
It is necessary to calculate the volume 



PdV= e~5- dV= e"«"(^j d€, 



"IlVJ''' 



dp^dpadqidq^dqs 



within the ellipsoid e = O, where for definiteness and simplicity the problem 
has been restricted to the case of the single system consisting of only three 
gongs which is the case first handled in the previous discussion. The easiest 
way of evaluating F is to make the change of variable 

Xi = \firaipi, Xj = \/va^Pi, x^ = s/vOsPs, x* = s/frv-^qi, x^ = \/v{qi - q{) 

I— v^ 

Ze = Vt(5'8 - qi), dpi- ■ ■ dqz = -jrfxj • • • dx^.f 

• As a Bystem with a definite length for the uppermost string is under discussion, 
the quantity v is supposed to be constant at this point; later it will be allowed to vary. 

t It may be recalled that to find the new element of volume when changing the variable 
in a multiple integral the Jacoblan is introduced. 
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The integral to be evaluated then takes the form 



1 rt,= c r 



where xl + xl + xl + xl + xl + xl^ C. 

This is a classic integral evaluable by means of the F— functions. The 
result is * 

^= 5 ^^^[^"'^^ = 1 C«r» = i ^v\ (7) 

Now the expression for the probability may be written 

PdV=e~^ . ^-,de; 
2 v^ 

and may be given a more convenient form by letting 

4> = ^og(-^) = 2 loge - 2 logv - log 2, 

so that I'dV= e^-^* de and T c^+* de = 1, 

the last condition being the new form of (2) • 

Now a second thermodynamic analogy may be set up. This time the 
analogy is between the quantities arising in the theory of a perfect gas and 
the average taken over a microcanonical ensemble. Equation (7) may be 
written as 

^v- = ke^oe y and compared with U^ V^ = Ke^. 

It is then seen that the following are analogical : 

In the Ensemble In the Gas 

the energy « the energy U 

the quantity v the volume V 

the quantity log V the entropy S 

* See Todhunter's Integral Calculus, p. 263, article 278. 
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The list may be further extended by the use of the fundamental thermody- 
namic formula 

d€ = dH-pdV= Td8-pdV 

compared with the result of differentiating e, namely, 

Then the following correspondences may be set up : 

In the Unsertible In the Gas 

(■^ tA = e-* V= ie the temperature T 

— I =- ) = ;r- = A the pressure p 

Equation ^v = 2c- * V Equation pV=2T 

Tliere seems to be no need to carry the table further. Similar results could 
have been found for the system with n gongs. 

This second analogy has the distinct advantage that the energy of the 
dynamical system corresponds to the energy of the gas, whereas in the earlier 
analogy it was the average energy over the whole ensemble which corresponded 
to the energy of the gas. This is not the place and the necessary fine dis- 
tinctions between closely related ideas in the theory of probability have not 
been sufficiently set forth to enter profitably into the discussion of the relative 
merits of the two analogies. On the whole, physicists seem to prefer the 
second, apparently mainly for the reason that the correspondence of the 
energies seems to them more natural ; but Gibbs evidently inclined rather to 
p efer the first and to feel that ( in view perhaps of the fects that for systems of 
so many degrees of freedom as a gas the energy of any one system could 
scarcely be expected to differ appreciably from the average energy and that 
we really do not know what the energy of the gas as a dynamical system is 
except in so far as it makes itself felt through the sort of average phenomena 
of pressure and temperature) the first correspondence was quite satisfactory ; 
there are, moreover, some other reasons why it is more satisfactory. 
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Conclusion. The foregoing simple dynamical system in which the com- 
putations of finding averages have been so easy as to be carried out without 
difficulty and for which the thermodynamic analogies are immediately compar- 
able with the ordinary theory of gases has been treated as far as possible in 
parallel with the general discussion of statistical mechanics given by J. 
Willard Gibbs in his Elementary Principles in Statistical Mechanics. That 
the simple and concrete case may perhaps be useful as an introduction to the 
general and abstract presentation of the statistical mechanics of any system which 
satisfies the Hamiltonian equations of motion is its only excuse for existence. 
There are, however, a number of matters which suggest themselves in connec- 
tion with the foregoing work and with statistical mechanics in general which 
may deserve a moment's notice. 

In the first place it might be considered remarkable that a system so 
simple and so apparently unlike a gas should give the thermodynamical equa- 
tions of a gas. A reference to Gibbs's work will show, however, that the 
fundamental differential equation of thermodynamics and the fundamental 
equation of statistical mechanics 

of C/ = Td8 - pd V and d7= @d(-~r))- Ad F or de = e-* Vd log V - Adv, 

whether in the canonical or microcanonical ensemble, will always con*espond 
when there is only one variable external coordinate in the system. It there- 
fore only remains to select the dynamical system in such a way as to give the 
proper special results for a gas, namely, the characteristic equation, the free 
energy, and so on. The selection was easy, especially in view of the general 
law of equipartition of kinetic energy. 

A few words may be said in explanation of. this famous law propounded 
by Maxwell and discussed by practically every student of the theory of gases 
since that time. In the case in hand, the kinetic energy is a sum 
7r{aipi + aipl + agpl) of the squares of the momenta and it is clear that the 
individual terms TrUipl, Tra^pl, ira^pl represent the kinetic energies of the in- 
dividual gongs. The average value of these partial kinetic energies is easily 
found. For instance 
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iraipi = 



£••/' 



Trample » dp^ ■ ■ • dq^ 
* — 'p—' 



/r •<• 'P 'n 
• • • / e « dpi- ■ -dq^ 

Traiple"^P'' dpi • • • e ® dp^ 



i — >ra,}>,2 — irOj P3« — «, 



rf^'s 



»o. Pi' — iraaPa' — «, 



e"hPi'dpi / . . . / e * (Zp2 . . . <Z?8 

Now the integrals though properly iterated, may be considered as merely 
multiplied : for^x does not occur in the last five provided the limits are not 
taken as functions of pi- Hence the last five integrals may be canceled from 
numerator and denominator and the calculation of the simple integrals remain- 
ing gives 

iruipl = J @. 

It therefore is seen that the average of each partial kinetic energy correspond- 
ing to each individual gong is ^ @ and the average value of the total kinetic 
energy is J 3© or Jn0. The fact that each gong has the same average kinetic 
energy as every other constitutes the law of equipartition of kinetic energy 
for this particular case. With the proper conventions the same result is ob- 
tained for any dynamical system. 

The proof that the average kinetic energy in a system of n degrees of 
freedom is ^ n 6 and that each degree of freedom has its aliquot share of the 
whole depends, as Gibbs gives it, on the assumption of the canonical distribu- 
tion including the hypothesis that the momenta may have any values 
between — oo and + so and that the limits of the integrations with respect to 
the coordinates qi, q^, •••»?»» ^o ^^ot depend on the momenta. If any of 
these assumptions were untrue, the proof would not hold. The assumptions 
seem ideally simple with the exception of the initial assumption of canonical 
distribution ; and at any rate the fact of the at least approximate truth of the 
law of equipartition of energy is so strikingly confirmed in the theory of gases 
as to justify the detailed discussion of any reasonable theory which has it as a 
consequence. This will perhaps serve as an answer to those who might inquire 
why the particular law of density or probability for the ensemble known as the 
canonical law was chosen for the work of this article ; there are moreover 
numerous more or less accurate and plausible "proofs" of the law or reasons 
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for it Avhich may be found in various treatises on the theory of gases and in 
Gibbs's book ; but they cannot be taken up here. Other assumptions for the 
density or probability would have enabled us to compute averages but would 
in all probability have been destitute of thermodynamic analogy. 

Although the general theory of statistical mechanics reached a high de- 
velopment at the hands of Gibbs, the detailed application of the method of 
investigation to systems that are of thermodynamic importance in physics has 
scarcely been touched ; special methods and approximations for obtaining 
averages, under the initial assumption not of the canonical distribution but of 
the law of equipartition of energy, have generally been employed. Indeed it 
is readily seen that to evaluate an iterated integral like 

/• • ■ dpy ■ ■ dp„j- ■ • e I -dqi • • • dq„ 

(even where the limits for the integrations with respect to the 5''s are taken 
as independent of the ^'s) is out of the question for most sets of limits and 
for most assumptions for the potential energy e, even if the simple assimaption 
be made that the kinetic energy consists of a smn 2a,jj| with constant coeiB- 
cients. For instance if Xi, yi, Zi and Xj, y^ Zj are the coordinates of the 
centers of two monatomic molecules of a gas, and if the usual law of attrac- 
tion inversely as the square of the distance and repulsion as the inverse fifth 
power of the distance be assumed, the potential energy would take the form 

'"^tr U(a;, - xjy + (y, - yjY + {z, - z^y 

+ (X, - x^y + (y* - yjY + (2.- - ZiYYl '• 

to which other terms must be added to keep the molecules within the vessel 
containing the gas. 

Despite the difficulties which seem to lie in the way of the systematic 
practical application of the statistical mechanics of a canonical distribution of 
dynamical systems, the fact that the theory leads to the fundamental differential 
relation of thermodynamics will probably lead to its detailed investigation. It 
is hoped that the presentation here given may be an aid to those who wish to 
make a start on this field of science. 

Massachusetts Institute op TECHHOi.oaY, 
Mabch, 1909. 



